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The rippling instability of a liquid sheet 
was first observed by Debr4geas, de Gennes, and 
B ro char d-Wy art [Science 279, 1704 (1998)] on 
a hemispherical bubble resting on a free sur- 
face. Unlike a soap bubble, it collapses under its 
own weight while bursting, and folds into a wavy 
structure which breaks the original axisymmetry. 
In fact, this effect occurs for both purely elastic 
and purely viscous (liquid) sheets, and an analogy 
can be made between the two mechanisms. We 
present a theory for the onset of the instability 
in both cases, in which the growth of the cor- 
rugation out of an inextensible initial condition 
is governed by the competition between gravita- 
tional and bending (shearing) forces. The insta- 
bility occurs for a range of densities, stiffnesses 
(viscosities), and sizes, a result which arises less 
from dynamics than from geometry, suggesting a 
wide validity. We further obtain a quantitative 
expression for the number of ripples. Finally, we 
present the results of experiments, which are con- 
sistent with our predictions. 

Every day, nature surprises us with structures 
and patterns of such beauty as to fill the scientist 
with wonder and the artist with envy. In the present 
paper, we address an instability which turns a uni- 
form, smooth, liquid bubble into a striking wrin- 
kled structure. The rippling of a deflating bubble 
was first observed by Debr^geas, de Gennes, and 
Brochard-Wyart [1]. In their experiment, 0.1 to 10 
cm 3 of air injected in a highly viscous liquid (r) ~ 10 3 
Pa*s) rises to the free surface, imprisoned in a hemi- 
spherical bubble of thickness e ~ 1 - 10 pm. If the 
latter is perforated at its apex by a needle, a circular 
opening expands exponentially at first. After about 
30 ms, the retraction slows down to a steady state. 
In the meantime, the air flow through the hole equi- 
librates the pressure difference, allowing the bubble 
to collapse under its own weight. As it deflates, an 
instability appears: the fluid sheet folds into a wavy 
structure, with radial ripples that break the original 
axisymmetry. In the absence of a detailed theory, 
Ref. [1] proposes a scaling estimate for the number 
of ripples n m [pgR 3 jK) , where p is the mass 

of the film per unit area, g the gravitational acceler- 
ation, R the radius of the hole, and K the effective 
bending rigidity of the sheet which was assumed to 
be elastic during the early stages of the rippling. Vis- 
cous flows eventually cause the corrugation to decay, 
as the hole retracts towards the thicker edge of the 


bubble. 

The rippling results from the competition be- 
tween compression, bending, and gravity. Each fluid 
element tends to fall under its own weight, but expe- 
riences a viscous resistance from its neighborhood. If 
the bubble were to collapse in a uniform, symmetric 
way, it would occupy a progressively reduced area, 
leading to an in-plane compression which would re- 
quire forces that far exceed the scale set by gravity. 
Instead, the film deforms in a nearly inextensional 
fashion by undergoing pure bending almost every- 
where. In terms of dynamics, this argument may be 
restated as follows; for a given (gravitational) force, 
the relative time scale associated with stretching is 
much larger than that for bending, and the surface 
therefore prefers to corrugate on short times, be- 
fore eventually relaxing into a uniform, thicker mem- 
brane. We note that this is valid only for creeping 
flows; the instability does not occur in the bursting 
of low viscosity films such as soap bubbles, owing to 
the rapidity of the retraction. 

This instability is reminiscent of buckling phe- 
nomena in slender elastic bodies [2]. In fact, creep- 
ing flows of liquid filaments and sheets may also lead 
to buckling (a striking everyday example being the 
coiling of a stream of honey when it reaches a toast 
[3,4]), an analogy we shall elucidate and exploit. For 
an elastic rod, buckling occurs at the longest possi- 
ble wavelength in order to minimize the bending en- 
ergy. In the bubble problem, however, gravity plays 
a distinctive role in determining the configuration. 
Though bending still favors large scale deformations, 
the gravitational energy is minimized for an almost 
flat sheet with as many tiny ripples as possible, and 
the optimal wavelength results from a compromise 
between the two. 

The rippling instability corresponds to the 
growth of the most favorable perturbation, trans- 
forming an axisymmetric initial condition into a cor- 
rugated pattern of a given wavelength. Although 
the bubble has the geometry of a sphere before col- 
lapsing, it is quite flattened by the time the ripples 
appear. For simplicity, we consider the unperturbed 
configuration to be a shallow cone, of slope a 1, 
described by its height above the surface, 

h = a(p 0 -p), (1) 
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where p is the cylindrical radial coordinate and p 0 
the radius of the base (Fig. 1(a)). As a further 
simplification, we do not incorporate the full burst- 
ing dynamics and corresponding viscous flows in our 
analysis, but carry it out with a specific hole size, 
which enters the theory as a parameter R. Once 
the hole grows steadily, the rate of retraction of the 
liquid film is given by the balance of viscous stress 
and surface tension (cr ss 20 mN/m), resulting in 
a constant velocity v ~ <7/77. It thus takes a time 
r ~ rjef <7 for the opening radius to increase by e. In 
this time, the liquid acquires a velocity V ~ gr due 
to gravity, larger than v by a factor V/v ~ 10 7 . Even 
if the liquid is viscoelastic, so that the retraction ve- 
locity is enhanced by a factor R/e (~ 10- 10 4 ) [5,1], 
the hole radius remains essentially constant while 
the instability occurs. 

To allow for the rippling of the surface, we 
perturb our truncated cone, without any loss of gen- 
erality, as 

h + 5h — a ( po — p) — 6a (p) 

+ £ S & n (/0 cos ( nd ) > ( 2 ) 

n> 1 


where 6 is the azimuthal angle. In the case of a 
thin elastic sheet loaded by its own weight, the two 
primary modes of deformation are in-plane stretch- 
ing and out-of-plane bending. If the scale of the 
deformation is of 0(1), the bending forces are pro- 
portional to Y e 3 /i, while the stretching forces are 
proportional to Yel [6]. Their ratio scales as ( efl ) 2 , 
so that for a given external force, here due to grav- 
ity, inextensional deformations are greatly preferred. 
For a thin, very viscous sheet, the primary modes of 
deformation are also bending and stretching (shear- 
ing) ones. Here however, the forces result from ve- 
locity gradients in the sheet. For deformations (of 
the bubble) of 0(1), the ratio of the viscous forces 
[7,8] due to bending, fae 3 // 2 ) dl/dt, to those asso- 
ciated with stretching, ( T]e)dl/dt , is (e//) 2 , as be- 
fore, largely favoring inextensional deformations if 
e « po/n m . (This condition is satisfied if the 
selected number of ripples n m is small compared to 
10 , which is the case as we shall see later) . Equiva- 
lently, for a given loading, the time scale correspond- 
ing to bending is smaller than that for stretching by 
a factor (e/l) , Thus, at the onset of the instability 
we may neglect perturbations that involve stretch- 
ing, both for the elastic and the viscous sheet. More 
specifically, we require (i) a vanishing Gauss curva- 
ture [9], leading to a developable surface, so that 

Sa (p) = Sa x (c - p ) , S(3 n (p) = 50 n x p (3) 

(where Sa , 50 n , and c are constants); (it) no in-plane 
strains [9], yielding 


53 (l + n2 ) (^") 2 = 4a5a. (4) 
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The perturbation thus consists of an axisymmetric 
flattening (—6a) along with ripples of magnitude 
S/3 n , analogous to the deformation of a tablecloth 
draping a circular table (Fig. 1(b)). We note that 
the inextensibility condition constrains the ripples to 
decrease in magnitude upon approaching the open- 
ing, in agreement with experimental observations. 
(The linear dependence of 6a and 5/3 n on p, though, 
remains to be checked ) 

To complete the characterization of the prob- 
lem, we have to specify boundary conditions for 5h. 
As the bubble rests on a liquid substrate, its edge 
cannot sink, so that 6h (p = po) = 0. Clearly, this 
cannot be satisfied by the simple ansatz of Eq. (3) , 
but requires one involving both bending and stretch- 
ing. However, this stretching is restricted to a small 
region near the edge [6,10], so that the inextensional 
assumption holds in the bulk, along with Eqs. (3) 
and (4). Although the condition 6h (p = p 0 ) = 0 is 
not compatible with Eq. (3), it is best approached 
by the choice c = p 0 , which also minimizes the grav- 
itational potential energy. 

In what follows we consider the elastic case 
first and then transpose our results to the viscous 
case. For inextensional elastic deformations, the en- 
ergy functional reads 


B[h + <5h] = / ^(surface) (gravitational energy 
J bubble 

+ bending energy) 

rPo 

= pdp de^i + Wh + shw 

Jr Jo 

{ fig(h + 6h) + ~[V 2 (h + 6h)] 2 }. (5) 


Here K = Ye 3 /12 (l — u 2 ) is the effective rigidity, 
where Y is the Young modulus and u the Poisson 
ratio. On substituting Eqs. (2-4) into Eq. (5) we 
obtain, to lowest order in the perturbation, 

5E = E[h + 6h] - E[h] 

= f^(f>£^") 2 {-7*V(f) 

n> 1 

+ [l- 7 i? 3 V'(^)]n 2 +n 4 } 

= ( 6 ) 

n>l 


where 



and 
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Here, we have used Eq. (4) to eliminate <5a, and 
the result is valid when (<5/ 3 n ) 2 ~ aJa « a « 1. 
7" 1 / 3 is an intrinsic length scale, arising from the 
competiton between gravity and bending elasticity. 

Each mode contributes an amount SE n to the 
change in energy, and rippling occurs if 5E n < 0 for 
some integer. In general, SE n < 0 for a range of dif- 
ferent n’s; the most negative variation corresponds 
to the maximally growing perturbation, and thus 
sets the wavelength of the instability. From Eqs. (6) 
and (7), we see that the instability is unexpect- 
edly suppressed when 7 < 7 C = j?~ 2 (po - iJm) _1 - 
Here, R m maximizes iiV(p 0 /i?), and is calculated 
from (p 0 - R m ) In (p 0 /Jl m ) = Rm^P (po/^m)* Fur- 
thermore, for a given bubble size po, rippling does 
not appear as long as R < Rc (p, K } p 0 ), where R c is 
defined by R\xp (po/J?c) = 7” 1 . Azimuthal continu- 
ity near a small hole translates into short wavelength 
deformations, associated with a high (bending) cost. 
Equivalently, if R < 7“^ 3 , one expects deforma- 
tions close to the hole to be forbiddingly expensive. 
Finally, the onset of rippling also depends on the ra- 
tio po/R (as seen from the implicit definition of Rc)- 
This is related to the fact that the wavelength of the 
deformation close to the hole is R/po times smaller 
than at the base, implying a bending cost ( po/R ) 2 
times larger, while the gain in gravitational poten- 
tial energy is proportional to (po — R). Minimizing 
SE in Eq. (6) yields the number of ripples as 



where Int{x} is the integer closest to x. This rela- 
tion improves on the estimate of Ref. [1], where the 
authors consider the short time elastic behavior, and 
establishes its domain of validity. 

The above calculation amounts to a linear sta- 
bility analysis involving the forces F = -5E/5h de- 
rived from the energy functional of Eq. (5). This 
balance of forces for the elastic sheet differs from 
the equations of motion of the viscous sheet only in 
that the elastic modulus Y is replaced by the viscos- 
ity 7} and the height h is replaced by the transverse 
velocity dh/dt of the sheet [8]. Indeed, for short 
times after the viscous bubble begins to collapse, this 
can be obtained by integrating the Stokes equations 
through the thickness [11,8]. Using such a proce- 
dure, one can write the viscous resistance to bending 


35 4( J(curv^ur e) ^ ^ ^ ^ effect j ve ben(J _ 

ing modulus of a liquid sheet is K = Ki = T]e 3 /3 t 
{y = 1/2 for an incompressible medium), where r 
is a time scale associated with the falling velocity. 
Thus, all the conclusions of the stability analysis for 
the elastic sheet, and in particular the expression 
for the number of ripples (Eq. (8)), may be trans- 
posed to the case of the bubble modulo a certain 
time scale related to the gravity-induced velocity of 
the fluid. Comparing the nascent ripples’ amplitude 
to the film thickness yields an estimate of this time 
scale as (e/p) 1 / 2 . However, as the bubble continues 
to fall, this time scale may vary and so may the num- 
ber of ripples. This last possibility is unlikely owing 
to the large forces required to introduce new ripples 
or remove existing ones. 

In order to check our results against experi- 
ment, we visualized the bursting of silicone oil bub- 
bles. Once the bubble is punctured with a sharp 
needle, its evolution is followed using a high-speed 
camera capable of recording up to 1000 frames per 
second. The resulting video is then analyzed to de- 
termine the radius po of the bubble, the hole size R 
at which the ripples are first observed, along with the 
number of ripples n*. As mentioned above, the hole 
expands very feist at first; by the time the pressure 
difference has equilibrated and the bubble begins to 
collapse, R is much larger than R c . In order to com- 
pare the experiments with the theory, in which R 
enters as a parameter, the latter is measured at the 
onset of the instability for each given size of the bub- 
ble. In Fig. 2, we plot n* as a function of po, which 
compares well with the theoretically predicted curve. 
Furthermore, the critical dependence of the rippling 
on the radius po is qualitatively confirmed by the 
experiments which showed a supression of the insta- 
bility for small bubbles. 

We conclude with a discussion of possible re- 
finements of the theory and their relation to the ge- 
ometric nature of the problem. A more complete 
theory would incorporate a (flattened) hemisphere 
as the initial condition, rather than a cone. Also, 
while we argued that the inevitable stretching is lim- 
ited to a negligible region, we have assumed that the 
latter is relegated to the lower edge of the bubble, 
disregarding the fact that the liquid is progressively 
drained out of the film by gravity so that the sheet 
is thinner around the hole and therefore more exten- 
sible (but also even more flexible) at the top than at 
the bottom. Though we expect, given the boundary 
conditions, that the stretching occurs quite generally 
near the base, improvements on the theory would 
be accompanied by the choice of more complicated 
ansatzes in Eqs. (1) and (2), and include a variable 
thickness e(p); the corresponding calculations will 
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be reported elsewhere [11]. Yet, the strong geomet- 
rical constraints involved in the problem are sugges- 
tive of the robustness of the results. The question we 
have answered is akin to that of applying a curved 
surface unto a flat one in the most economical way, 
a problem which has taxed cartographers for many 
centuries and lies at the birth of differential geom- 
etry. It is also somewhat of an inverse counterpart 
to the problem of fitting a flat sheet to a three di- 
mensional landscape, which has been studied in var- 
ious contexts [12-15], and is an issue that still vexes 
fashion designers. The relevance of the geometrical 
constraints is manifest, for example, in the strong 
dependence of the rippling on the size of the open- 
ing, which is closely related to a well-known theorem 
due to Gauss [16], Jellett [17], and others, accord- 
ing to which (loosely put) a closed surface cannot be 
bent without being stretched, while an open surface 
can be bent inextensionally. Similarly, we find that 
a smaller hole implies a relatively stiffer bubble, and 
hampers the rippling. While the detailed form of the 
functions ip and xp arise from the physical constraints 
and dynamics imposed by the the energy functional 
(or forces) and the various boundary conditions, the 
essence is in the geometry. 
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FIG. 1. Stroboscopic images of a collapsing liquid 
bubble of size is po = 1mm and thickness e ss 100/xm. 
The silicone oil has viscosity t) = 10 3 Pa s, surface tension 
cr = 21mN/m, and mass density 0.98g/cm 3 . (a) 30ms af- 
ter the film is punctured by a sharp needle, the bubble 
shows a retracting hole of radius R = 1.4mm at its apex, 
but no ripples yet. (b) 30ms later, the bubble loses its 
axisymmetric shape. The radius of the hole remains es- 
sentially constant, at R = 1.6mm, while the ripples grow. 
(The decrease in amplitude upon approaching the hole 
is noticeable.) The inset displays a schematic side view' 
of the essentially conical deflating bubble at the onset of 
the instability, with the important quantities involved in 
the phenomenon. The extreme shallowness allows for a 
perturbative treatment in the slope a of the cone. 
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FIG. 2. Plot of the number of ripples n* as a func- 
tion of the bubble radius po, comparing the experimental 
measures (points) with the theoretical prediction (solid 
line). These data were gathered using solicone oil of vis- 
cosity T) = 600Pa * s and bubbles of thickness e % 30pm. 
The errors in the measurement of po arise from menis- 
cus effects which are more important in smaller bubbles. 
The bursting time elapsed up to rippling is measured to 
be of order of 1 to 5 times (e/g) 1/2 , consistent with our 
proposed mechanism for the formation of the corruga- 
tion. For each experimental realization, the ratio po/R 
was measured at the onset of the instability, and the 
corresponding dependence of R on po was used to get 
a theoretical curve n m = n’ (po). The latter is plotted 
here for a slope a « 3 0 (ss O.OSrad) of the cone, consis- 
tent with our perturbative treatment, and in agreement 
with direct observation. The dashed line represents the 
best fit of the scaling form n* ^ (pgR z /K) ^ [1], where 
R is chosen as the relevant length scale. If R is replaced 
by po, the above expression for n* may be closely fitted 
(up to an overall multiplicative factor) to our predicted 
curve, showing that the size of the bubble is the domi- 
nant length scale within the present experimental range 
and conditions. We also note that the predictions com- 
pare best with the experimental data for larger bubbles 
(e/po < 1), for which an inextensible model is better 
suited. The inset shows a top view of the fully devel- 
oped ripples, from which n * is measured. 
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